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ABSTRACT 



As a foundational element describing relativistic reacting waves of relevance to 
astrophysical phenomena, the Rankine-Hugoniot relations classifying the various 
propagation modes of detonation and deflagration are analyzed in the relativistic 
regime, with the results properly degenerating to the non-relativistic and highly- 
relativistic limits. The existence of negative-pressure downstream flows is noted 
for relativistic shocks, which could be of interest in the understanding of the 
nature of dark energy. Entropy analysis for relativistic shock waves are also 
performed for relativistic fluids with different equations of state (EoS), denoting 
the existence of rarefaction shocks in fluids with adiabatic index T < 1 in their 
EoS. The analysis further shows that weak detonations and strong deflagrations, 
which are rare phenomena in terrestrial environments, are expected to exist more 
commonly in astrophysical systems because of the various endothermic reactions 
present therein. Additional topics of relevance to astrophysical phenomena are 
also discussed. 

Subject headings: hydrodynamics: shock waves — ISM: kinematics and dynamics 



-3 - 



Introduction 



The classic al theory o 



fluid dynamics ( jWilliams 



combustion waves is well established in the study of reactive 



1985 



Zel'dovich et al. 



1985 



Law 



20061 ). Different from 



hydrodynamic shock waves which have been extensively studied, supersonic detonations 
and subs onic deflagrations are sustained by reactio ns in the fronts of fluid discontinuities 



(see, e.g. 



Landau fc Lifshitz 



1959; 



Williams 



19851 ) . and as such are expected to yield rich 



varieties of fluid dynamical responses. 



Axford and Newman ( lAxford 



1961 



Newman fc Axford 



19681 ) first adopted the 



concept of reactive fluid dynamics in ast rophysics by considering the 



hydrogen ionization fronts around stars. 



Blandford fc Ostriker 



(1978 



dynam ics of the 



1980) advanced 



the mechanism of particle acceleration through astrophysical shocks such as those in 
supernova systems, which in principle should also bear the dynamics of reactive flows. 



Indeed, it is now commonly believed that nuc 



support the explosion of 



2000 



Gamezo et al. 



2003 



vpe. 



ear def 



agration and detonation waves 



1969 



Hillebrandt fc Niemeyer 



a supernovae (lArnett 
20051 ). recognizing nevertheless that the transition mechanism 
from deflagration to detonation is still not clear. Additional studies of astrophysical 
processes that have suggested/invoked combustion processes include relativisti c detonation 



waves as a possible mechanism for the false vacuum decay (ISteinhardt 



for the microwave background fluc tuations ([Gibson 



fireball model for the 7-ray burst (IMeszaros 



19821 ). relevant 



20051 ) in the early universe; the 



20021 ) involving relativistic 7 photons and 



electron/positron pairs; and the unique role of the ISM shocks in vari ous chemical process es 



in mo 


ecular clouds, particular r 


f in early-phase 
) as observed in 


2007 




Flower & des Forets 


2010 


Furthermore, a recent study ( 


Gao & Law 


2011) 



McKee fc Ostriker 



201 ll ) on the evolution of supernovae remnants 



has also incorporated reaction to account for the accelerative expansion of the Crab Nebula. 
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It is thus clear that integration of combustion theory offers rich potential in the study of 
astrophysical phenomena. 

To adopt reactive fluid dynamics in the study of astrophysical systems, a general form 
of com bustion wave theory in relativistic fluids is needed. While theories of relativistic sho ck 



waves (ITaub 



19481 ) and detonation waves in highly- relativistic fluidqj (Steinhardt 



idJ3 



have been advanced, and the rel ativistic shock waves for radiating fluids (jCissoko 



and for fluids in magnetic fields (IMallick 



1982) 



1993) 



20 111 ) have been presented, a general analysis of 



combustion waves describing subsonic deflagration as well as supersonic detonation waves, 



for all relativistic fluids, has not been performed. Consequent 



we shall int egrate the essential featu res of relativistic fluids (ILandau fc Lifshitz 



Anile 



19891 ) and combustion waves (IWilliams 



1985 



Law 



y, as a first, necessary st ep, 



1959; 



20061 ). within the context of 



the Rankine-Hugoniot relations, and identify the various possible relativistic combustion 
waves and their properties. The relativistic theory of deflagration and detonation waves 
has also been studied in quark-gluon plasmas using the bag equation of state for the 



quark matter (IGyulassy et al. 



1984 ), which is a useful reference for the present study of 



19571 ) with the adiabatic index T for more general 



combustion waves in a Synge gas (ISynge 
astrophysical applications. 

Since the relativistic shock is an important component in relativistic flows, it will be 
separately studied first in Section 2, in which an interesting solution involving negative 
pressure in the downstream fluid is identified, and entropy analysis to the shocks are 
founded. In Section 3, the general relativistic form of deflagration and detonation waves is 
presented, with its proper degeneracy to the non-relativistic and highly-relativistic limits. 
Detonation and deflagration waves, as well as special cases of Chapman- Jouguet waves and 



1 Here highly-relativistic fluids refer to those hydrodynamic systems whose speed of fluid 
particles are very close, or equal, to the speed of light. 
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isobaric waves are analyzed therein. The possible existence of weak detonations and strong 
deflagrations in relativistic astrophysical environments is then discussed, which is followed 
by summary of the present work, in Section 4. 



2. Relativistic shock waves 



We first briefly outline the fundamentals of relativistic fluid 



of relativistic shock waves (ITaub 



1948 



Landau &; Lifshitz 



1959 



dynamics an d the theory 



Liang 



1973). 



2.1. Basic equations governing the relativistic gas 



Based on the energy-momentum tensor in the local rest frame of a relativistic fluid 

^ e ^ 



p 
p 
y p J 

momentum and energy conservations across a shock front are: 



(1) 



Uiul +pi = u 2 u\ + p 2 



(2) 



and 



Wl7lMl = W 2 72^2, 



(3) 



respectively (ILandau fc Lifshitz 



1959 



Steinhardt 



19821 ) . Here e is the fluid energy density, 



which includes the rest-frame energy nm^c 2 of the fluid particles and the specific internal 
energy e (see Equation flB])), p is the pressure of the fluid, u> = e +p the specific enthalpy 
per unit volume, mo the rest mass of one particle and n the particle number density 
which can be the number density of any charge (e.g., baryon number or lepton number) 
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that is conserved across the shock front. The four-velocity u of the fluid has the form 
u = (3j, where (3 = v/c is the velocity in unit of the speed of light and 7 = 1/(1 — /3 2 ) 1 / 2 
is the Lorentz factor. Subscripts 1 and 2 denote upstream and downstream quantities, 
respectively. The particle number conservation 



niUi = n 2 u 2 



(4) 



is a complementary condition of continuity across the shock. 

For analysis pertinent to the stellar interior, the ISM or intergalactic media (IGM), 



and the early universe, the relativistic ideal gas is described b y the fo 
EoS) with the adiabatic index Y (i.e., the Synge gas, cf. 



state 



1982 



Cissoko 



1993): 



Synge 



lowing equation of 



1957 



Lanza et al. 



r = i+*, 

where p = nrrtQ is the rest-frame mass density and e the specific internal energy (cf. 
Equation (jSJ)). Then the fluid energy density can be written as 



(5) 



e = p(c 2 + e) = nrriQC 2 + 



1 



r- 1 



P, 



(6) 



which can also be considered as the definition of the specific internal energy e. Using the 
EoS (jSJ), the specific enthalpy has the following form (IKennel fc Coroniti Ill984j ). 



u = e + p = nm c + 



(7) 



The EoS with 4/3 < T < 5/3 (ITaub 



1948 



Anile 



19891 ) is a general form for relativistic 



fluids. In the highly- relativistic regime (r = 4/3), the internal energy greatly exceeds the 
rest-frame energy of the particle, i.e., e ^> c 2 , and t h e flui d is radiation dominant with 



V 



\pe = ~e. This is the EoS used in ISteinhardt 



(1198a ! in the study of the bubble 



growth during the false vacuum decay in the early evolution phase of the universe. In 
the non-relativistic extreme (r = 5/3, cf. Anile 1989), e c 2 , the EoS assumes the form 
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p = |pe. Comparing this expression with the classical ideal gas law p = pkT, with k 



being the Boltzmann constant, the internal energy is just the kinetic energy e = |/cT for a 
monatomic gas in non-relativistic fluids. Another special EoS with T = 2/3 accounting for 
the dark energy in the universe will be discussed in Section 2.3. 



The sound speed in relativistic fluids is given by (cf. 

'dp 



Landau fc Lifshitz 



19591 ) 



v s = c\ 



de' 



(8) 



whose four- velocity is u s = /3 S 7 S = — (1 — (^) 2 ) 1//2 . Following the above discussion, we have 
v s — cf \/3 in the highly-relativistic regime and v B = 



I 2 in the non-relativistic regime. 



2.2. The relativistic shock adiabat 



By defining a variable x = uj/n 2 for fluids in both sides of the wave front, we readily 
obtain the expressions for the particle flux (j — ri\U\ — n 2 U2) by considering momentum 
conservation ([2]) and particle flux conservation (j3J) 

-f= P -^, (9) 

X2 — X\ 

jy additionally conside ring the energy conservation ([3]) (cf. 



and the so-called shock adiabat 



Taub 



1948 : 



Landau &: Lifshitz 



1959 : 



Steinhardt 



19821 ): 



x 2 uj 2 -xiui = (p2 -pi)(x 2 + Xi). (10) 

Referring to the form of the specific enthalpy (J7|) under the EoS, the variable x can be 
expressed as 



x 



— (nm c + 



TP 



2 . .— y 

+ 



mlp 



'IV. 



u' \ r — v ) p r — i p 2 

Introducing Equation (TTTj) to the particle flux dHJ), we obtain the Rayleigh relation for 
relativistic fluids, accounting for the conservation of number density and momentum: 



-it 



c 2 (V-l) 



-{pV 2 - 1) 



(12) 
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Here V — 1/p is the specific volume, and notations with hats are reduced variables, i.e., 



V = P2/P1, V = V2/V1 an d c = c/ y/piVi. In the limit of c — > 00 and u-\ —> 0, the 



Landau fc Lifshitz 



Jlayleigh 



19591 ): 



relation degenerates to that in the non-relativistic regime (cf. 

Ml = pi-, (13) 

r(v-i) 

where M — vj sjY'pj p is the non-relativistic Mach number; while in the limit of c — > 0, it 
assumes the highly-relativistic form: 

(r-i)(p-i) 



"1 



T{pV 2 - 1) 



(14) 



By introducing x given by ( Ulj) to the shock adiabat IjTOj) , we derive the Hugoniot 
relation for relativistic shocks: 



r + i 
r-i 



(pV-l)-(p-V) 



r 



(fv* - 1). 



(15) 



r-r v_ ' ' (r-1) 2 

The Hugoniot relation additionally takes into account of energy conservation ([3]) and as 
such describes all possible discontinuities across a shock by considering all conservation 
equations. In the non-rela tivistic extreme of c — > 00 and U\ — > 0, Equation ( fl5|) reduces to 



(^Landau fc Lifshitz 



19591) 



p + 



1 



r - 1 



4r 



(16) 



r + uv r + iJ (r + 1) 2 ' 

In the highly-relativistic extreme of c — > 0, the RHS of Equation ( fTBT) vanishes, yielding 
The highly-relativistic shock adiabat, i.e., Rayleigh and Hugoniot relations, have been 



derived and discussed in, e.g., 



Steinhardt 



(119821 ). 



Hugoniot lines represented by Equation ffl5l) are illustrated in Figure [T] for non- 
relativistic, relativistic and highly-relativistic fluids. Along reference line a, it is seen that 
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the density increase is higher for relativistic fluids than for non-relativistic fluids for the 
same pressure enhancement ratio across the compression shock. And along reference line 
b the pressure reduction is smaller in relativistic fluids than in non-relativistic fluids for 
the same density dilution ratio across the rarefaction shock. These differences are due to 
the fact that pressure assumes a larger proportion of the total energy in relativistic fluids 
when the flow speed gets closer to the speed of light, such as the downstream of rarefaction 
shocks along reference line b. 

Typical shock solutions for relativistic and highly-relativistic fluids are shown in 
Figure [2J One distinct difference from non-relativistic shocks is that the Rayleigh lines 
for relativistic fluids are not straight lines any more, which can be easily seen from the 
form of Equation (TP2]) . A more essential difference between relativistic and non-relativistic 
shocks is that compression shocks (p 2 > Pi, V% < V\) can only be achieved with very high 
upstream flow speeds in relativistic fluids as indicated for example by the state {u\ = 1, 
v i = -^c) in Figure |2j Furthermore, when the upstream flow speed is reduced by half 
{ui = 1/2, vi = -^gc), only rarefaction shock solution (p 2 < pi, V 2 > Vi) exists. The 
criterion distinguishing these two types of shocks is the sound speed of the upstream flow, 
i.e., u\ = A= and v\ = 4=c, when the only tangency state between the Rayleigh and 
Hugoniot lines is the (1,1) point. This classification of shocks and the criterion are the same 
as those for the non-relativistic shocks. 



2.3. Negative pressure downstream flows 

Under certain conditions of particle flux (|9]), the downstream pressure can assume 
negative values. Take the highly-relativistic case as an example. Combining equations ( I14p 
and (fTTl) . and by setting the adiabatic index to V — 4/3, we readily obtain the relation 
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between the reduced downstream pressure and the upstream four-velocity: 

ti? = i(3p+l). (18) 

Then for upstream velocity in the range < U\ < the reduced downstream pressure is 
negative, i.e., — | < p < 0. Furthermore, under this condition the reduced specific volume 
is an imaginary number according to the highly-relativistic Hugoniot relation (117j) . i.e., 
V 2 < 0. That is, for shock waves in highly-relativistic fluids with a normal upstream flow 
(pi > 0, Vi > 0) whose speed U\ is smaller than a negative pressure state can be 
achieved in the downstream flow with the specific volume (or density) being an imaginary 
number. The p — V 2 diagram showing the highly-relativistic shocks with negative-pressure 
downstream flows (u\ = 0.2 and 0.3) is shown in Figure [3j The physical interpretation of 
this type of negative pressure fluids with imaginary density (defined as Type I) is not clear 
so far, although its existence can be ruled out based on entropy considerations, as will be 
demonstrated in the sequel. 

Another type of negative pressure fluids (Type II) is the case of T = 2/3 in the EoS 
©, i.e., 

P=~\p^ ( 19 ) 
which can be interpreted on the basis of dark energy in the universe. Under this equation 
of state, the combination of highly-relativistic Rayleigh ( FH| and Hugoniot ( TTTj) relations 
readily leads the form of the upstream four-velocity: 

«?=i(-30+l). (20) 

For the upstream flow speed u\ > the reduced pressure is negative, accounting for a 
negative pressure downstream fluids with the upstream fluids being normal. This type of 
negative pressure fluids has a real-number specific volume (or density) according to equation 
f|T7|) . Shock waves involving Type II negative pressure fluids are shown in Figure HJ 
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Entropy analysis showing the availability of these two types of negative pressure fluids 
are presented in the following. 



2.4. Existence of rarefaction shocks 



Alls 



rock waves should 



relation ( iLandau fc Lifshitz 



ollow the law of entropy increase. Applying the thermodynamic 



1959|) 



d 



n 



Td( a - 



n 



dp 



(21) 



in the weak shock wavqj, we have the following form through Taylor expansion: 



n 2 U\ \n 2 n\) ri\ 



Here T is the temperature, a the entropy per unit proper volume and s = a/n the entropy 
per particle. From equation (TlO|) . we have another expression of ^ — i.e. 



CU 2 OJl 1(1 1 . 
= T, — + — (P2-P1J- 

n 2 ni 2 Vn 2 ni - 



(23) 



It should be noticed that in achieving equation (f23|) . the higher order term — (p 2 — Pi) 
has been omitted. Combining equations (122]) and (123]) . we readily obtain the relation of the 
entropy difference across the shock front: 

Vn 2 ni/ 2 Vn 2 ni/ 2 V / s 6 V c 

By introducing the expansion of — with respect to p 2 — Pi 

1 1 fd(l/ri)\ , , l(d 2 (l/n)- 



{p 2 - Pl f. (24) 



n 2 ni 



( d(l/n) \ 

\-^p-) s ^-^ ' 2 



dpi 



(P2 - Pi) 



(25) 



2 The weak shock wave assumes the discontinuity in every quantity to be small. 
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to equation f[24|) . we have the simplified expression of the entropy difference across a 
relativistic shock: 

^_!!l = S2 _ Sl = JL(^))( P2 _ Pl )3. (26) 

The law of entropy increase requires that s 2 — s\ > 0. 

From the equation of state (J5j), the particle number density can be expressed as 

- = (L^IW. (27) 
n p 

By using this expression in the entropy equation (l26l) . we have 

a 2 cTi (r-l)m ei 

= — — (p 2 - Pi), (28) 

from which it is easily identified that for normal upstream fluids with p\ > and e\ > 0, 
rarefaction waves with p 2 — p\ < only exist for the adiabatic index T < 1. So under the 
regime of weak shock, the Type I negative pressure fluids with T = 4/3 does not exist 
due to the violation of entropy increase law; while Type II negative pressure fluids with 
T = 2/3 is expected to be a real physical existence. However, since it is still questionable 
as whether the entropy analysis in weak shocks can be generalized to normal strong shocks, 
the existence or non-existence of Type I negative pressure fluids in strong shocks needs to 
be further studied. 



3. Relativistic detonation and deflagration waves 

When exothermic or endothermic reactions are involved in fluids, detonation or 
deflagration waves can be excited, which renders the dynamics different from those of 
non-reactive flows, as will be analyzed in the following. 
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3.1. Rankine-Hugoniot relations for relativistic combustion waves 

Representing the overall energy release per unit mass by q, which is positive and 
negative for exothermic and endothermic reactions respectively, the specific enthalpies per 



unit volume across the reaction front 



are 



u)i = e\ + pi + pq and 

lo 2 = e 2 + p 2 (29) 

for upstream and downstream flows, respectively. Following the procedure in deriving 
Equations ffT2l) and f|T5|) . the Rayleigh and Hugoniot relations for relativistic reactive flows 
are given by: 



(p-l) = -u{ c\V - 1) + ^(PV 2 - 1) - q , (30) 

r _ ^ r t + v 

)q^ 

(31) 



c 2 



r + i 



T (pV-l)-(p-V)-2q] =^p(l-V 2 )-^ f - I ^(p 2 V 2 -l) + (p+ f ± I )g + g 2 



where the reduced heat release is q = —q. 



pi 



Representative Hugoniot lines described by relation (|3T|) are shown in Figure |5j We 
see that while the non-reactive Hugoniot line passes through the (1,1) point in the p — V 
diagram, the Hugoniot lines for flows with exothermic and endothermic reactions are above 
and below the non-reactive one, respectively. This difference leads to the possible existence 
of weak detonations in endothermic reactive flows, which will be discussed further when 
considering Figured 

We next note that Rayleigh lines in relativistic reactive fluids do not pass through the 
(1,1) point in the p — V diagram (see e.g. Figure [6]) as they do in non-reactive relativistic 



3 The reaction front is normally slim in dimension relative to the entire flow under con- 
sideration. 
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flows (Figure [2]). This is because in relativistic reactive flows, the reaction heat release 
not only is part of the total energy, but it is also present in the momentum conservation 
equation (see Equations ([2]) and (|29p ). which then leads to the presence of the q term in the 
Rayleigh relation (!30l) . 

We now separately discuss the solutions for the relativistic detonation and deflagration 
waves. 



3.2. Detonation waves 

Detonation wave solutions for relativistic fluids with exothermic and endothermic 
reactions are shown in Figure E] and Figure [3, respectively. 

Figure [6] shows that detonation solution exists for exothermic reactive fluids with 
relatively large upstream speeds, i.e., u\ > \/2; while there is no detonation solution for 
fluids with relatively slow upstream speed, i.e., u\ < l/y/2. However, the criterion of the 
upstream flow speed for the existence of detonation is higher than the speed of sounc3, 
which is different from the case of non-relativistic fluids, for which this criteria is exactly the 
speed of sound. This is due to the existence of the heat release q in the Rayleigh relation 
fl30|) for relativistic reactive flows. By numerically exploring Rayleigh lines with different 
upstream speed ui, we find that the criteria for the existence of detonation wave solutions 
are u\ = 1.0 and u\ = 0.9 for relativistic (c = 1, T = 3/2) and highly-relativistic (c = 0, 
T = 4/3) gases, respectively, both with the energy release q = 1 (see Figure |6]). Another 
interesting feature is that, while flows with u\ = a/2 have both strong and weak detonation 
solutions, higher speed flows with u\ = 4 have only weak detonation solutions. The reason 



4 In the highly-relativistic fluids, the speed of sound is u s — 1/ v2; while in normal rela- 
tivistic fluids, this value is lower. 
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for the disappearance of strong detonation for exothermic fluids with high speed upstream 
flows is that the reaction heat q is greatly amplified by the large value of the upstream flow 
speed Mi (this is a relativistic effect), which then leads to large pressure compression ratio p 
to satisfy the energy conservation requirement involved in the Hugoniot relation (l3Tj) . 

Detonation solutions in endothermic reactive fluids (Figure [7]) are quite different. 
Only weak detonation can be found, and its solution exists for various upstream flow 
speeds ranging from subsonic to supersonic. The reason for the extensive existence of 
weak detonations is that endothermic fluids tend to absorb the fluid kinetic energy of 
any strength in order to initiate the endothermic reactions (e.g., ionization of the ISM), 
which then behave as detonations with relatively lower pressure and density compression 
ratios. In Figure [7J the Rayleigh lines for all upstream flows with different speeds have 
single intersections with the Hugoniot line, which shows that there is no threshold for the 
existence of weak detonations. On the other hand, strong detonations cannot be formed for 
flows with any upstream speed. 



3.3. Deflagration waves 

Deflagration wave solutions for relativistic fluids with exothermic and endothermic 
reactions are shown in Figure |8] and Figure [U respectively. 

Figure E] shows that exothermic reactive fluids have different kinds of deflagration waves 
for different upstream flow speeds. Specifically, when the upstream flow speed is higher 
than U\ = 0.5, there is no deflagration solution in both relativistic and highly-relativistic 
flows. Taking into consideration that no detonation solution exists for exothermic flows 
with upstream speeds U\ < 1/a/2 (see Section 3.2 and Figure [6]), we conclude that for 
a certain range of the upstream flow speed (e.g., 1/2 < u\ < l/\/2 in Figure |6] and [8]) 
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in relativistic exothermic reactive fluids, neither detonation nor deflagration waves can 
form. More specifically, numerical criteria for the existence of deflagration waves have been 
identified as U\ = 0.37 and u\ = 0.39 for the relativistic and highly-relativistic gases in 
Figure [HI with the energy release q = 1. Figure [8] further shows that, for the relatively 
low upstream speed of U\ = 0.3, both strong and weak deflagrations exit; but when the 
upstream speed becomes much lower, i.e., Ui = 0.2, strong deflagration cannot be achieved, 
leaving only the weak deflagration. This can be understood by referring to the Rayleigh 
relation fl30|) . which shows that as the upstream speed u\ decreases, the pressure rarefaction 
ratio p (< 1) should be higher for the same reduced specific volume V. Then at a certain 
low value of u\, there is no intersection between the Rayleigh and Hugoniot lines for higher 
V, i.e., strong deflagration does not exist. 

The wave response is however quite different for endothermic reactive fluids, as shown 
in Figure |9j No deflagration solution can be found for flows with upstream speeds of 
u\ = 0.2 and 0.5. Deflagration waves emerge only when the upstream speed is as high 
as u\ = 1.0 (higher than the sound speed). The numerical criteria for the existence of 
deflagration waves are u\ = 0.7 and u\ = 0.6 for relativistic and highly-relativistic gases 
in Figure [9j respectively, with the energy release q = — 1 . Referring to the deflagration 
waves with low upstream speed in exothermic reactive fluids, we expect that endothermic 
deflagrations need higher upstream flow speeds to propagate. Furthermore, even when the 
upstream flow speed is higher than the speed of sound, shock cannot form because part of 
the fluid kinetic energy is absorbed by the endothermic reaction, which is exactly the case 
of Ui = 1.0 in Figure M It is noted (see Section 3.2 and Figure [7]) that the existence of weak 
detonations for a large range of the upstream speed is a distinct feature of endothermic 
flows. Consequently weak detonation waves, instead of deflagrations, should be the expected 
dynamic structure in endothermic reactive flows. 
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3.4. Chapman-Jouguet waves 

As a parallel analysis to the non-relativistic Chapman-Jouguet waves, the tangency 
point of the Rayleigh and Hugoniot lines, representing a limit case of the wave solutions, is 
analyzed for relativistic fluids. At the tangency point the slopes of Rayleigh and Hugoniot 
lines are equal to each other, so slopes of both curves in the p — v diagram are calculated 
here according to Equations (!30|) and (13T1) : 



fdp\ c 2 {p-l) + 2^ l pV{p-l) 

\(IV) Rayleigh & (V ~ 1) + (V 2 ~ 1) - q ' 

/#\ + 1) + h^pV_ + 2 w ± W p 2 V 

W^Hugoniot " C 2 (^V - 1) - ^(1 - V 2 ) + 2 J ^ w pV 2 

And the solution to the relation 

fdp\ = / dp\ 

\(IV' Rayleigh V dV ' Hugoniot 



(32) 



(33) 



(34) 



shows the criterion of the waves represented by the tangency point. Typical numerical 
solutions for the Chapman-Jouguet waves as criteria for the existence of detonation waves 
are given in Section 3.2 and Figure [6] captions. 

It is noted that in the non-relativistic extreme of c — > oo, the tangency solution 
degenerates to the normal Chapman-Jouguet wave, i.e., 

p - l r^iP + 1 



v-i 



(35) 



Law 



20061 ) . which leads to 



This condition implies that the downstream flow is sonic (c.f. 
the classical Chapman-Jouguet wave in which the downstream flow does not affect the 
the wave front so that the wave propagates steadily. However, no simple solution of the 
tangency point can be found for relativistic fluids as the heat release q are involved in both 
relations (1321) and (BSD. 
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3.5. Isobar ic waves 

A comparison between Figures [6] - [9] and Figure [2] shows that Rayleigh lines for 
exothermic and endothermic fluids intersect at different points from those of non-reactive 
fluids in the p — V diagram. Specifically, in exothermic reactive fluids, the Rayleigh lines 
intersect at a point with p = 1 and V > 1; while for endothermic reactive fluids, the 
intersection is at p = 1 and V < 1. This difference between reactive and non- reactive fluids 
is caused by the decrease and increase of flow densities in constant-pressure exothermic 
and endothermic waves, respectively. The example of Figure [7] shows that because the 
intersection of Rayleigh lines for endothermic reactive fluids is to the right of the Hugoniot 
line, one can easily achieve weak detonation solutions. 

Since the intersections of Rayleigh lines for reactive fluids with different reaction heats 
q have the same reduced pressure p = 1, we expect the existence of a common solution 
with p = 1 for any given reaction heat q. Then the Rayleigh relation (I3(jp can be readily 
converted to the following form by considering (p — 1), instead of p, as a common factor: 

-u\ \^V 2 + c 2 V - (c 2 + q + ^) 

(p - 1) = l - s — j -. (36) 

1 + ujV 2 ^ V ' 

If we calculate the (positive) V root of the numerator in the RHS of this equation, the 
coordinate of the intersection of the Rayleigh line is 

T - 1 



p = l, V 



2T 



-c 2 + Jc 4 + 4-J—[c 2 + q + 



1/J 



(37) 



In the non-relativistic limit, c — > oo, the intersection is at (p = 1, V = 1); while in the 



highly-relativistic limit, c — > 0, the intersection is at (p = 1,V = J 1 + ^jr-q). (See Figures 
|6]- |9]for the variation of the intersection.) In the limiting case of small reaction heat in 
highly-relativistic fluids, i.e., q = ^q 1, the V axis of the intersection can be expressed 
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as V = 1 + ^2F"<?, which can be readily converted to the form of 

r-io 

V 2 -Vi = — (38) 

The above equation implies that in an isobaric wave, the specific volume increases and 
decreases for exothermic and endothermic reactions respectively. In other words, the flow 
density decreases in an exothermic reaction as the fluid releases energy and increases in an 
endothermic reaction due to the absorption of energy, if the pressure is kept constant. 



3.6. Further considerations of weak detonations and strong deflagrations 



Based on entropy considerations, raref action shocks norma. 



the a diabatic index T > 1 (cf. Section 2.4, 



Landau fc Lifshitz 



y do not exist for fluids with 



1959 



Zel'dovich fc Raizer 



19661 ). so those weak detonation as well as strong deflagration waves involving rarefaction 
shocks also normally do not exist. However, in non-relativistic fluids, weak detonation 
and strong deflagration waves can be found in systems with en dothermic rea c tions or 



1961 



Williams 



phase transitions, in which a rarefaction shock is not necessary ( lAxford 
19851 ). While these understandings still apply in exothermic relativistic detonations and 
deflagrations, i.e., the weak detonation solutions in Figure [6] and strong deflagration 
solutions in Figure |8] should normally not exist in realistic systems, there is an exception 
in that for very high upstream flow speeds (e.g., Rayleigh lines with u\ = 4 in Figure |6]), 
weak detonation waves could exist as no rarefaction shock structure is required in such 
detonations. Instead, these exothermic weak detonations with high upstream speeds are 
similar to the weak detonation in endothermic reactive fluids (Figure [7]), as both reactions 
can be directly ignited in the shockless waves because the kinetic energies of the upstream 
flow is large enough to initiate the reaction. It is therefore reasonable to expect that 
weak detonation as well as strong deflagration waves of this kind are more common in 



high-speed astrophysical flows, and they are expected to have different structures from the 
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tra ditional Ze 
(cf. IWilliams 



dovi ch - von Neumann - Doring (ZND) structure in strong detonation waves 



19851 ). One such example is that of relativistic endothermic reactive fluids, 
for which weak detonation waves instead of weak deflagration waves are the dominant 
dynamics as can be seen from Figure [7] and Figure [H 

We further note that the non-relativistic form of endothermic combustion waves 



has been studied for interstellar gas ionization ( lAxford 



196ll ). The ionization of neutral 



hydrogen atoms is an endothermic reaction which forms an 'ionization front' in the 



interstellar media. This front, separating the unionized (H I) and the fully ionizec 



(h in 



regions, is dynamically identical to the detonation (or deflagration) front. lAxford I (Il96ll ) 
has shown the existence of both weak detonation and strong deflagration waves theoretically 
and what we discussed here is a generalization of these results to the relativistic regime. 
Besides the ionization of neutral hydrogen atoms around hot stars, the r e-ionization (also 



of ne utral hydrogen atoms by stars and quasars) process of the universe (IMiralda-Escude 



20031 ) can also be considered as an endothermic combustion wave. Another endothermic 
fluid process in astrophysics is the inverse Compton sc attering, the mechanism for a large 



varie ty of high energy X-ray and 7-ray radiations (e.g., 
201lk 



De Jager fc Harding 



1992 



Tavani 



For relativistic fluids with T < 1, rarefaction shocks are allowed through entropy 
consideration (cf. Section 2.4), and as such lead to the formation of weak detonations. This 
is obviously another way of having weak detonations in astrophysical relativistic fluids. 



4. Conclusions 

As a theoretical foundation to study the dynamics of astrophysical systems, the basis 
of relativistic reactive fluid dynamics is formulated here in terms of the Rankine-Hugoniot 
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relations for detonation and deflagration waves in normal relativistic reactive flows. The 
relativistic shock theory is revisited and the mathematical solutions of two types of 
negative pressure downstream flows are identified. Normal relativistic and highly-relativistic 
detonation and deflagration wave solutions are constructed for both exothermic and 
endothermic reactive flows. The existence of endothermic reactions in astrophysical 
phenomena extends the family of potentially realizable combustion waves in terms of weak 
detonations and strong deflagrations, which do not exist extensively in terrestrial situations. 
These theoretical results can be applied in astrophysical systems such as supernovae 
explosions, 7-ray bursts and the false vacuum decay in the early universe, for which the 
conventional shock theories do not hold. 
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Fig. 1.— Hugoniot lines for non-relativistic (Equation ffl6|) with T = 5/3), relativistic 
(Equation (fT5l) with c = 1 and T = 3/2) and highly-relativistic fluids (Equation ( TTT1) with 
T = 4/3). The ordinate is the downstream pressure relative to the upstream one; the abscissa 
is the downstream specific volume relative to the upstream value. All curves pass through the 
point (1,1). We notice that for the same downstream to upstream pressure compression ratio 
(reference line a), the density increase is higher for relativistic fluids than for non-relativistic 
fluids; while for the same downstream to upstream density diluent ratio (reference line b), 
the pressure reduces less in relativistic fluids than in non-relativistic fluids. 
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Fig. 2. — Rankine-Hugoniot shock solutions for relativistic (c = 1, T — 3/2) and highly- 
relativistic (c = 0, T = 4/3) fluids in the p — V diagra m. The Rayleigh lines are not straight 



lines as they are in the non-relativistic case (see, e.g-. lLandau fc Lifshitz II1959I ). We can see 
that for (highly-)relativistic fluids, compression shocks with P2/P1 > 1 can only be achieved 
for relatively higher upstream fluid speeds (here u\ = 1 for example). And for relatively 
lower upstream fluid speeds (here U\ = 1/2 for example), only "rarefaction shocks" with 
P2/P1 < 1 (which normally do not happen due to entropy increase) exist. In between these 
two kinds of shocks is the sound-speed-upstream-flow, i.e., the Rayleigh line with U\ = l/y/2, 
which does not have any shock solution. 
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Fig. 3. — Shock solutions for highly-relativistic (c = 0, T = 4/3) fluids in the p — V 2 diagram 
with negative pressure downstream flows. The square of the specific volume is also a negative 
value, which means that the downstream specific volume (or density) is an imaginary number 
if we assume the upstream flow to be normal fluids. We can see that there is no intersection 
between the Hugoniot line and the Rayleigh line with u\ = 0.5, while intersections exist for 
Rayleigh lines with u\ = 0.2 and u\ = 0.3. In fact, for highly-relativistic fluids, only when 
the upstream flow speed is in the range of < u\ < ^ (see Section 2.3) can we achieve 
negative pressure flows via this kind of shocks. However all this shocks do not follow the 
entropy increase law (cf. Section 2.4) so this Type I negative pressure fluids physically do 
not exist. 
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Shocks in highly— relativistic fluids with 
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Fig. 4. — Shock solutions for highly- relativistic (c = 0, T = 2/3) fluids in the p — V 2 
diagram with negative pressure downstream flows. Shock solutions exist for Rayleigh lines 
with upstream speeds in the range of u\ > i.e., u\ = 0.5, u\ = 1.0 and u\ = These 
Type II negative pressure fluids follow the entropy increase law (cf. Section 2.4) thus are 
physically available. 
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Fig. 5. — Hugoniot lines for relativistic gas (c = 1, T = 3/2) and highly-relativist ic gas 
(c = 0, T = 4/3). In the first quadrant, the Hugoniot line with positive reaction energy 
(here q — 1 for example) is above the shock Hugoniot line with q = 0, and the Hugoniot line 
with endothermic reaction (here q — — 1 for example) is below the shock Hugoniot line. It 
is also noted that Hugoniot lines are also present in the forth quadrant, hence allowing the 
existence of shock or reactive wave solutions with negative pressures fluids. One such shock 
solution has been demonstrated in Figure HI 
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Fig. 6.— Detonation wave solutions for relativistic gas (c = 1, T = 3/2) and highly- 
relativistic gas (c = 0, T = 4/3) with an exothermic reaction (q = 1). For relatively 
low-speed fluids with ui = l/\/2, there is no intersection between the Rayleigh lines and 
the Hugoniot lines with q — 1, i.e., no exothermic detonation wave exist. For intermediate- 
speed fluids with ui = y/2, binary intersections can be found between the Rayleigh lines 
and the Hugoniot lines, implying the existence of both strong and weak detonations for 
the exothermic reactive flows. For high-speed fluids with ui = 4, there are only single 
intersections between the Rayleigh lines and Hugoniot lines, which are weak exothermic 
detonations. The criteria for the existence of detonation waves are u± = 1.0 and u± = 0.9 
for the left and right panels, respectively, as obtained from numerical explorations. It is also 
noticed that for both the left and right panels, all Rayleigh lines intersect at a point with 
q = 1 and V > 1. 
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Fig. 7.— Detonation wave solutions for relativistic gas (c = 1, T = 3/2) and highly- 
relativistic gas (c = 0, T = 4/3) with an endothermic reaction (g = —1). The ordinates are 
in the logarithm coordinates. For fluids with all upstream speeds shown here (ui — 1/ v^2, 
u\ = \/2 and ui = 4), there are single intersections between the Rayleigh lines and Hugoniot 
lines with q = — 1 in the detonation region (p > 1), which means that weak endothermic 
detonations always exist. In both panels, Rayleigh lines intersect at a point with q — 1 and 
V < 1. 
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Fig. 8. — Deflagration wave solutions for relativistic gas (c = 1, T = 3/2) and highly- 
relativistic gas (c = and T = 4/3) with an exothermic reaction (q = 1). For relatively 
high-speed fluids with u± = 0.5, there is no intersections between the Rayleigh lines and the 
Hugoniot lines with q — 1, i.e., there is no deflagration waves. For intermediate-speed fluids 
with ui = 0.3, binary intersections can be found between the Rayleigh lines and Hugoniot 
lines with q — 1, implying the existence of both strong and weak deflagrations for the exother- 
mic reactive flows. For low-speed fluids with u\ = 0.2, there are only single intersections 
between the Rayleigh lines and Hugoniot lines with q — 1, representing strong exothermic 
deflagrations. The criteria for the existence of deflagrations in the upper and lower panels 
are U\ = 0.37 and U\ = 0.39, respectively, as obtained from numerical explorations. In both 
panels, Rayleigh lines intersect at a point with q = 1 and V > 1. 
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Fig. 9. — Deflagration wave solutions for relativistic gas (c = 1, T = 3/2) and highly- 
relativistic gas (c = 0, T = 4/3) with an endothermic reaction (q = —1). For low-speed flows 
with ui = 0.2 and u± = 0.5, there is no intersection between the Rayleigh lines and Hugoniot 
lines with q = — 1 in the deflagration region (p < 1), which means that endothermic defla- 
grations do not exist. For relatively-high-speed flows with u± = 1.0, there are intersections 
between the Rayleigh lines and Hugoniot lines with q = — 1 in the deflagration region (p < 1), 
representing the existence of endothermic deflagrations. Numerical explorations show that 
the criteria for the existence of deflagrations are U\ = 0.7 and Ui = 0.6 for the upper and 
lower panels, respectively. In both panels, Rayleigh lines intersect at a point with q = 1 and 
V < 1. 



